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Abstract 

The form factors for semi-leptonic B decays, B — ► irlvi, are calculated under collinear fac- 
torization approach. The end-point divergences are regularized by a £-regularization, where £ 
means the collinear fraction of the spectator anti-quark of the B meson. The form factors 
are calculated up-to 0{a s /mB)- The complete Oil/ tub) contributions from the tt meson are 
calculated explicitly by a collinear expansion method. A well-defined power expansion scheme 
is built such that the Ofl/ms) contributions are about 30% of the leading order contribu- 
tions. A small value i*+(0) = 0.164 is found. This confirms the SCET result i<+(0) = 0.17 
from B — ► -kit decays. The form factors are calculated for < q 2 < 16GeV 2 , where q 2 is 
the invariant mass of the lepton pair in B — > irWi. An extrapolation of the form factors to 
q 2 > 16GeV 2 is made to obtain \V ub \- 2 J 26A2Gey2 dq 2 (dT / dq 2 ) th = (5.71 ± O.^ps' 1 . We 
determine \V ub \ = (3.95 ± 0.13 oxp ± 0.32 th ) x 10~ 3 from the world averaged branching ratio 
Br(B -> irl&i) = (1.36 ± 0.09) x 10~ 4 . 
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I. INTRODUCTION 



B factories have obtained precise results on Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements of the standard model (SM) [1]. The LHCb has already started running. More re- 
sults on CKM matrix elements with higher accuracy are expected in near future 2J. Being an 
important element of CKM matrix elements, \V u b\ still contains large uncertainties dominated 
by theoretical ones. For example, the \ V u b\ from inclusive B —>■ X u Wi processes contains 10% 
uncertainties, where the 7% uncertainty comes from the 60 MeV uncertainty for m;, and the 



other 3% uncertainty comes from experiments. ThejK 
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contains a 10 — 15% uncertainty (j 

from the B — > 7r form factor, F? n , and a 6% uncertainty from experiments [20j, 
addition, the q 2 spectrum of B — > txIv\ has been well constrained experimentally [22( 

To take full advantage of the experimental precision for exclusive B — > ttIvi, it is necessary 
to pin down the theoretical uncertainties on to few percentage level. However, it is 
still a difficult task. Currently, nonperturbative methods, includin g Q CD light-cone sum 



rule (LCSR) ,Q and lattice QCD (LQCD) Q, GJQ, Q 



are available. 



Innovative parameterization (PA 
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methods with model independent inputs from theories 
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Naive application of perturbative 



or linear end-point divergences 



for the end-point divergences 
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(pQCD) to Fr* needs to account for logarithmic 
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. Different methods have been proposed 
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361,0. In principle, the 



pQCD method can give a model independent determination for the |V^&| and its precision 
can be improved order by order by perturbation theories. Based on factorization theorem 



38j, the pQCD expression for the form factor can be written into a factorized formula in 



terms of hard scattering function and nonperturbative meson distribution amplitudes (DAs). 
The hard scattering function contains short distance contributions and can be calculated, 
perturbatively. The meson DAs contain long distance contributions and are universal. Once 
the meson DAs are determined from other processes, the factorization formula can make 
model independent predictions. However, due to large uncertainties associated with the 



pQCD form factors 
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351 ]. the pQCD method has not been applied to derive 



equally precise |T4&| as the other methods, such as, LCSR, LQCD, and PA. In this paper, 
we would like to improve the precision order of the pQCD form factors, such that one can 
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derive the \V u b\ with equally theoretical uncertainties as the other approaches. 

There are two compatible pQCD methods, the collinear factorization denoted by C* , and 
the kx factorization, or , PQ CD factorization, denoted by k T . The k T has been widely used 
for F BlT 



31 



32 



33. 



34 . l3q . In k T , end-point singularities are supposedly solved by the 
intervention of parton transverse momenta 39 1. However, the transverse parton momenta 

n 

would induce large logarithms a s \n kx from higher loop corrections [311 ] . In addition, there 
are also large logarithms a s In 2 x associated with subleading twist (twist-3) contributions 



3l| . These large logarithms need to be re-summed to be Sudakov factors [3l|. The Sudakov 



factors are expected to suppress the contributions from end-point re gions [331] . In practical 



applications, some criteria for the Sudakov factors are needed 



31 



33, 



391 ] . When k T is 

generalized to include subleading order contributions in the I/tub expansion, the subleading 
order, 0(1 /uib), corrections dominate over the leading order contributions 3l|, |33j. Intrinsic 
transverse degrees of freedom of the meson wave functions (for B meson and pion ) are needed 
to cure the ill behavior of the power expansion 351 ]. 

Unlike the complicate features and related issues of k T , C$ has a simple structure and 
is directly related to the parton model (PM)[4Q]. It is expected that, once the end-point 
singularity can be regularized within C* , the C* formalism would be instructive for both 
heory and phenomenology. In the approach proposed by Akhoury, Sterman, and Yao (ASY) 



26l | , the heavy quark effective theory (HQET) and Sudakov re-summation were incorporated 



with C* . The end-point divergent problem is solved in the ASY approach. However, a 
dynamical zero point was found and a small partial decay rate for B — > -nlv\ was obtained. 

In this paper, we propose a different approach to solve the end-point divergent problem 
and avoid the problems in the ASY approach. The key solution is a £-regularization (denoted 
by £ R ) which can effectively regularize the end-point divergences. The £ R has been applied 
to effectively regularize the end-point divergences in twist-3 hard spectator and annihilation 
contributions in charmless hadronic B decays 411] • In this paper, we show that £ R is also 
effective for F Rn . The suppression of end-point radiative corrections is provided by the B 
meson distribution amplitude (DA). This provides a stronger suppression effect than any 
Sudakov factors. The extension of application range of q 2 is given by including subleading 
corrections in 1/wig expansion. The complete 0(1/ rug) contributions from the 7r meson 



side are calculated by using a co 



which is developed by Yeh 41 



42 



linear expansion (CE) method for exclusive processes, 



431 ] . The linear end-point divergences in the 0(1/ tub) 
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contributions are solved by a simultaneous use of the £ R and a non-constant twist-3 pseudo- 
scalar (PS) DA. The factorization of the B — > n form factors has been shown valid under 
{3l| and soft collinear effective theory (SCET) 36], respectively. There lacked explicit 



regularization methods for the end point divergences in these previous proofs [31|, [36j. Our 
approach provides practical calculations for the form factors to show the factorization up to 
0{a s /mB)- The non-constant twist-3 PS DA exists for a 7r meson in its energetic state. On 



the other hand, a constant PS DA is usually used in the literature. As shown in j4l|, the 
constant PS DA is appropriate to describe a 7r meson in its chiral state (, or, a soft pion), 
but not an energetic pion. Unexpected large contributions associated with the constant 

33] and LCSR Q, Q, H, 0, 44] . If the constant PS 



31 



PS DA are already noticed in 
DA corresponds to the soft pion state, then these large contributions from the constant 
pion PS DA can be identified as soft dominated contributions instead of hard dominated 
contributions as expected in the calculations performed in the k T and LCSR. This point of 
view of using the pion state corresponding to the DA to identify the kinds of contributions 
(soft or hard) is different from the traditional way of using the scaling of the relevant 
contributions. This provides another method to identify the considered contributions. The 
unexpected large contributions from the pion PS DA in k^ T and LCSR could be overestimated. 
Since end-point singularities at leading and subleading order in l/mg expansion can be 
effectively solved, the F R q are calculable under . Another effect associated with £ R is 
that, for u ~ 0(1), the divergences from rj ~ 0(A/i71b) can be regularized. This extends the 
application range of G* from small q 2 ~ to moderate q 2 ~ 16 GeV 2 . Note that the relevant 



energy scale of a s is set as t = 1.65yl — q 2 /m 2 B , which is about 1 GeV at q 2 = 16 GeV 2 
and the coupling constant a s is about 0.48, or, a s /ir ~ 0.16 < 1. There are also possible 
subleading order contributions in l/ms expansion from the B meson side. Only I /tub 
contributions from the two parton Fock state of the B meson are calculated. In summary, 
we plan to make the following progresses in theory. 

1. The is applicable for F R $ at leading twist order by £ R . 

2. 0(1 /tub) two parton contributions are shown calculable under C' . 

3. The twist-3 three parton corrections from the pion are explicitly calculated. This is a 
first result in the literature. 
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4. The complete twist-3 contributions are shown less than the twist-2 contributions for 
< q 2 < 16 GeV 2 . A well-defined power expansion is given. 

The result is applied to extract \V u bF Rn (0)\ and \ V u b\ from the world averaged branching 
fraction Br^ for semi-leptonic decays B —>■ -nlv\. A fitting method is used to determine 
the value of a parameter uob for the B meson distribution amplitude. A parameterization 
form, Ff(q 2 ), of F Rn (q 2 ) versus q 2 is given by a minimal x 2 fitting. Our analysis gives 
\V u b\ = (3.95±0.13 exp ±0.32 th ) x 10~ 3 with experimental and theoretical errors. This agrees 



well with the world averaged value of \V u b\, \V u b\ = (3.95 ± 0.35) x 10 3 [451] . The fit form 
factor F c (q 2 ) predicts F? n (Q) = 0.16 which confirms the founding F? CET (0j = 0.17 by the 



soft collinear effective theory (SCET) from an analysis for B — > txtx decays 36]. 

The organization is as follows. (, R is defined and shown effective for leading twist contri- 
butions in Section II. The comparison between the l; R and the fc-r-regularization (denoted 
by k R in this paper) is given in this Section. To generalize C* for higher twist contributions, 
the collinear expansion method is used in Section III. In Section IV the F R q are explicitly 
calculated up-to 0(a s /mb). In Section V, the form factors are applied to determine \V u b\ 
from experiments. Last two Sections are devoted for comparisons and discussions. Two 
Appendices are given. 

II. LEADING TWIST B -> vr FORM FACTORS, END-POINT DIVERGENCES, 
AND ^-REGULARIZATION 

In this section, we first review how the end-point divergent problem of the leading twist 
B — > 7r transition form factors can arise. We then define the £ R and explain how it is 
effective for end-point singularity. The B — > it form factors F Rl ^(q 2 ) are defined by 

(n( P7T )\qYb\B(P B )) = 2FfVK + l^tf) - (F^(q 2 ) - F ^(g 2 )) ^ ~ , (1) 

where q = Pb — Pn- Another set of form factors, fi£, is also used in literature. Their 
definitions are 

(*(p)\<rfb\B(P B )) = fi n (q 2 )P^ + fi n (q 2 )^ ■ (2) 
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F+q and are related by the following identities, 

F^(q 2 ) = \{f^{q 2 ) + ft{q 2 )), (3) 
F^{q 2 ) = \((2 - v)fx m (q 2 ) + r]ft(q 2 )) , (4) 

where r\ = 1 — q 2 /m? B . Under q 2 — > 0, F^{q 2 ) = F^ n (q 2 ). At maximal recoil limit (the 
energetic limit for the 7r meson), the form factors F+ n (q 2 ) and P B7r (g 2 ) becomes identical to 
cancel the q 2 — > pole in Eq.([T]). The B meson's momentum P B is defined in the B meson's 
rest frame P B = (tub, 0,0,0) = (P B ~,P B ~,Q±) with P B = P B = ms/V^. The n meson 
momentum p n is written as p^ = \{j]m B , 0, 0, i]m B ) = (p+,0,0j_) where p+ = i]m B /\/2. 
q = Pb ~ Pw = (m B (l - T)/2), 0,0, -r)m B /2) = (q + ,q~,0±) where q + = (1 - T])m B /V2 
and q~ = m B j\pl. The light-cone coordinate system will be used in this work. Under 
C% partons carry collinear fractions of external meson momenta. The spectator anti-quark 
of the B meson carries a momentum 1>L = (0,l~ p ,0±) where lj p = £P B = ^m B /\^2. The 
b quark's momentum is defined as Pjf = (P 6 + , P 6 ~, 0±) with P b ~ = £Pg~ = ^m B j\[2 and 
P fc + = m 2 /(v2£m B ). The b quark is defined on-shell. This is different from the usual 
treatment in the literature that the b quark is assumed almost on-shell, P 2 ~ m 2 B . The quark 
q inside the w meson is defined to carry a momentum l q = (/+, - , 0j_) with Z+ = ur]m B /y/2. 
The anti-quark q inside the 7r meson has a momentum lq = (7^~, 0~, 0j_) with lq = ur]m B /\/2. 
£ and u are momentum fraction variables and £ = 1 — £ and u—l—u. E = m B j \[2 is used 
in the following text. 

Under large recoil condition, 77 — »■ 1, the 7r meson has an energetic momentum p^ ~ 
m B /y/2 ^> m n . We assume that the virtual gluon carries a hard-collinear energy scale. 
The PQCD is applicable because the involved strong coupling constant a s (t) at the inter- 
action energy scale t = y/i]A h m B with A h = 0.5 GeV is around 0.3. Under C?, the twist-2 
contribution to the matrix element (ttI^^I-B) is written as 

M tw2 ^ = fwfB——^—— f d&TiO f dutf{u)H tw2 ^ u) , (5) 
iV c Jo J 

where H tw2,tl (^,u) = P~( a )' tw2 ' M (£, u) + i^ 6 )' tu,2,,1 (£, it) denote the hard scattering functions 
for the lowest order Feynman diagrams depicted in Fig. [1(a) and (b). 0^ 2 (£) and </%(u) 
are the B meson's and it meson's leading twist LCDA, respectively. f B and f n are the B 
meson's and 7r meson's decay constants. Cp = (N 2 — 1)/ (2N C ) and N c are color factors. a s 
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is the strong coupling constant. H^ a ' b ^' tw2 ' fl (^,u) are written as 



U 



(6) 



U 



t]£u(£ — rju)E 



9 fir 



(7) 



One can arrive at the form factors fi^' 1 ™ 2 



(<f) 



Tra s (t) Cf 
E 2 W c 

na s (t) Cf 
E 2 ~N f 



■Mb J d&Tit) J dutf(u) 



(1 + r])r]u - (1 + r])£ + r] 

T] 2 £u(£ - T)U) 



(8) 
(9) 



where rj — 1 — rj and u — 1 — u. The form factors F^' tw2 from fi2' twA by Eqs.(j3ED are 
written as 



[■Bir,tw2 



v — £ + v u 



7 Bn,tw2 



2E 2 ~N, 



Ub J d&Ttt) I du^{u) 



r) + (1 - 2r))£ - (1 - 2r])r]u 
r)£u(t; - T)u) 



(10) 



(.11) 



(j) n (u) is linear in u. /j 8 ^'*™ 2 are finite. However, if (£ — r/u) is approximated to be (—rju), 
then / 2 B7r '* 1 " 2 becomes logarithmic divergence at u — > 0. This is the end-point divergent 
problem of the B — > 7r transition form factors. The key point is that, without loss of 
generality, the denominator of the internal b quark propagator is approximately to be, (Pb + 
k) 2 — ml ~ rn 2 B (C, — i]u) + 0(A QCD A) + 0(Aq CD ) where A = m B — m b , £ ~ 0(A/m B ) and 
A = Aqcd- The error terms are of next-next-to- leading order in 1/wig expansion. In the 
end-point region of u, u ~ 0(Aqcd/"1b), £ is as large as rju for r] ~ 0(1). £ retained in 
(Pb + k) 2 — ml is necessary. 



Akhoury, Sterman, and Yao (ASY) 26] proposed an approach to combine the heavy quark 
effective theory (HQET), a resummation of Sudakov double logarithms, and C* for exclusive 
processes. In the ASY approach, the b quark of the B meson is almost static at the energy 
scale much less than the b quark's mass, mj. The exchanged gluon in the Feynman diagram 
in Fig. [1] (b) carries soft energy. The b quark inside the B meson can be effectively described 
by an effective field, h v , with v = P B /m B of HQET. According to HQET, the light degrees of 
freedom of the B meson are identified as a brown muck such that only soft interactions can 



W (6) 

Figure 1: One gluon exchanged Feynman diagrams for the hard scattering functions of the C* 
amplitude for the B — > ttWi decays. The external meson states \B) and (ir\ are not shown. The 
cross vertex denotes the vector operator 7^. The diagram (a) describes the one gluon exchanged 
process between the bottom quark b and the spectator anti-quark q of the B meson. The diagram 
(b) describes the one gluon exchanged process between the quark q and the anti-quark q of the ir 
meson. 

exist between the spectator anti-quark and the b quark inside the B meson. There needs a 
subtraction operation to separate long distance and short distance contributions from the B 
meson side. The subtraction operation is equivalent to neglect £ in the denominator and £ in 
the numerator of the H^' tw2 '^ in Eq.(J7|). It is obvious that once £ and £ are subtracted from 
the H^' tw2,fl term, the end-point singularity in /rf^'*™ 2 is solved. When loop corrections are 
concerned, there would arise large infrared (IR) logarithms from the internal b quark line 
in the Feynman diagram in Fig. 1(b). A resummation over the large IR logarithms gives 
Sudakov form factors. Although the ASY approach can successfully solve the end-point 
divergent problem, there exist two problems. The first one is that there is a dynamical zero 
point at 77 = 1/2 in the \fi + / 2 |. We find that the dynamical zero point is due to the neglect 
of the £ factor in the H^' tw2 '^ term. Refer to Eq.([7j). The £ factor is recovered as one 
include the dynamical contributions from the b quark. This can be understood by referring 
to Appendix |A] The dynamical zero point is avoided in £ fl . The second is that small partial 
rates for B — ► ttWi were predicted in the ASY approach. This implies that the ignored 
contributions in the ASY approach are important. In £ R , the ignored contributions in the 
ASY approach are recovered. See Section V for this point. In our approach, two problems 
in the ASY approach are solved. 

In kJp, the divergent term (rju) -1 is replaced by (k 2 -. + rjuE 2 )' 1 . While the transverse 
momenta of partons are remained in the parton propagators, the partons become off-shell 
because only collinear and transverse momenta are kept. Off-shell partons may radiate 
infinite soft gluons as they pass through space before they compose into external mesons. 
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Re-summation of soft gluon radiations results in Sudakov factors, which are expected to 
suppress the soft radiations from the end-point region. However, some criteria are required. 
This is an uncertainty of 

£ is of order 0(A/E) that the B meson distribution amplitude <Pb(0 has a peak at 
£ ~ 0(A/E). In the end-point region of u, u ~ 0(A/E), <p^(u) ~ 0(A/E). The overall scale 
of the leading part of the form factor f^ n ^ tw2 i s counted as 0(1/ E 2 ) times the decay constants 
fnfB- Similarly, y^-*™ 2 j s dominated by large u ~ 0(1) and is counted as 0(1/ E 2 ) times 
the decay constants f^fs- Although f^ n ' tw2 anc l /vf 71 '*™ 2 receive contributions from different 
configuration regions, they are of the same order. Accordingly, contributions of subleading 
order in I/tub expansion should be included. They are calculated in next Section. 

The 7r meson's twist-2 DA <p^(u) is modeled to be its asymptotic form <f)^(u) = 6uu by 
neglecting its scale dependence. This is fulfilled for the precision of 0(a s ). The B meson's 
twist-2 DA 0^ 2 (£) is assumed to be modeled as 

<Pb (0 = -^- ex p( (I 2 ) 

(jJ B Ll>b 

which satisfies the following conditions 

d^B 2(I \0 = 1, (13) 



o 



v ^B m (0 ™b 



The value of A# satisfies the condition 6A# < 4A with A = mg- mj,. Different models for 
the (J) 1 ^ 2 (£) have only tiny differences if they are constrained by Eq. ffT^l) [47 ]. We note that 
the integration range over £ in the calculation of the form factors is [0, 1] instead of [0, oo]. 
The difference between these two integration ranges is equal to 



poo pi 

A = / d^(0 - / d^(0 
Jo Jo 

m B -. \ r m Bl 

+ 1 exp 

Ld B J u B 

~ 1.29 x 10 ( ~ 4) , 

where the last estimated number is calculated by using u>b = 0.46 GeV and rriB = 5.28 GeV. 
A is negligible. 

To distinguish from the for end-point divergences, we name the retain of £ in the 
internal b quark propagator as the £-regularization, It is instructive to see how £ R works 



for end-point divergences. Let's consider the 77 — > 1 limit of F^' tw2 {q 2 ). Except of the 
relevant parameters, the most singular part is the integrations 



Jo Jo £«(£ - 



U) 

It is seen that once £ is neglected in (£ — u)" 1 , an end-point divergence arise to be log(w). 
The retain of £ results in, for the integration over u, 



I dutf(u) 1 = 6(l + £ln|-z7r£). 



One can observe that the original logarithmic divergence logw is indeed regularized by £, 
effectively. The further integration over £ is analytic because the end-points for £ — > or 
£ — > are prevented by the .B meson distribution amplitude 0^ 2 (£)- F° r example, the 
integration for the most singular term is analytical as 

/Ws 2(/) (0flnf = -[l-e^ + (a-l)( 7s + r(a) + lna)] , 
Jo s s 

where a = tub/ujb, 1e is the Euler number, and T(a) is the Gamma function. Another 
widely used model for 0^° 2 (£) has the form 



" ^ K 2 + fB (i-?) 2 ) 2 



:tw2(II) 



where iVg and e# are determined by Eq. ffTBl . Although the integration can not be expressed 
explicitly, the result is also analytic because the 0^f 2 ' 7/ has a stronger suppression effect on 
the end-points, £, £ — > 0. Of course, the final result would depend on the model for <^ 2 (£). 
This model dependence is due to our rare knowledge for the B meson. However, these two 
models give almost the same numerical results for p^' tw2 j n For example, we can use 
ujb = 0.46 GeV to obtain 

F+ ir,tu ' 2 ' (2) (0) = 0.110, 

F Bn,tw2,(II)^ =0.104, (14) 

where ^^ 7r > tu ' 2, ( / )(o) and p , ^> tw2 '( 11 ) (q) are calculated by using 0^(£) and 0£^(£), respec- 
tively. iVs = 0.1536 and e# = 0.0061 are used in the calculation for j7^ 7r >*™ 2 >v r - r )(o). 

This shows explicitly that the £ i? indeed regularizes the end-point divergence logw. The 
suppression of contributions from the end-point region in £ fi is provided by the B meson's 
DA, 0^° 2 (£). The success of £ fl implies that the dynamical role of the spectator (anti-)quark 
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of the B meson is important in the studies of heavy to light processes. This is contrary 
to the heavy to heavy processes, such as B — > Dlu, in which heavy quark symmetry is 
useful. In the heavy quark infinite limit of the heavy to heavy processes, the light degrees 
of freedom of the heavy mesons remained independent of the h v — > h v /lv process, described 
by a Isgure-Wise function, ((v ■ v'). h v ( v >) are effective fields for the b and c quarks of the B 
meson and the D meson, respectively. If the same picture applies for B — > irlv processes, 
then the b — > ulv transition should be in a similar condition as that of h v —>■ h v >lv process. 
This means that the final state ir meson should be in its soft pion state, at which the pQCD 
method is inapplicable in principle. In fact, the light degrees of freedom of the B meson 
would experience violate fluctuations during the b — > ulv transition proceeds. This implies 
the applicability of pQCD [48j]. The heavy quark symmetry is inapplicable in B — > nlvi 
decays. We will see later that the intervention of the £ variable becomes a straightforward 
step under CE. See the next Section. 

III. COLLINEAR EXPANSION 

In this section, we describe how CE can be applied to derive twist-3 two parton and 
twist-3 three parton contributions. In this work, we consider only the twist-3 three parton 
contributions from the 7r meson. The three parton contributions from the B meson are 
complicate and left to other places. We begin with the amplitude for the diagrams in Fig. [2J, 

M " = /^y»/(^ Tr|i/ " (!B ' !) * B(iB) *' (i)l (15) 

+ ■■■ , 

where the first line corresponds to the four parton interactions (for the diagram depicted in 
Fig. [2(a)), and the second line represent five parton interactions (for the diagrams depicted in 
Fig.[2Jb)). The higher Fock state interactions are neglected. The H^^Ib, I), and H^(Ib, I, I') 
describe the parton interactions in the hard scattering center. The $s(Zs), $ir(l), $^(/,/'), 
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(a) (b) 



Figure 2: Feynman diagrams for the hard scattering functions and H„ of the C' amplitudes for 
the B — > ixlvi decays. The external meson states \B) and {tt\ are shown by bubbles with symbols 
B and 7r, respectively. The central bubble represents parton virtual interactions between relevant 
partons from the B meson and pion. The cross vertex denotes the vector operator 7^. The diagram 
(a) describes the interaction processes with four partons involved. The diagram (b) describes the 
interaction processes with five partons involved. There are the two quark partons and one gluon 
parton of the pion taking participate in the interactions. 

are the hadron functions denned as 

®b(Ib) = J d 4 ze aB Z (0\q(z)b(0)\B), (16) 

<M0 = J rf 4 ^(7r|g(^)g(0)|0), (17) 

$£(M0 = f d^ze^e^ 1 '- 1 ^' {ir\q(z)(-g)A v (z')q(0)\0} } (18) 

where color and spin indices are not shown explicitly. The gauge links between particle fields 
in the matrix elements are implicitly understood. The trace operation Tr is taken over the 
spin and color indices. 

The factorization of the amplitude into partonic and hadronic parts is composed of three 
steps, the factorization of loop parton momenta, the color index factorization, and spin index 
factorization. The factorization of loop parton momenta is performed by means of a Taylor 
expansion for the partonic part and followed by relevant integral transformations. The color 
index factorization and the spin index factorization are similar and can be done by color 
algebra and Fierz identities, respectively. These three steps are shown explicitly below. 

Once the partonic part is separated from the hadronic part of the amplitude, it is an 
important task to examine whether the partonic part suffers from soft divergences. This is 
the proof of the factorization theorem. The 0(a s ) analysis for the twist- 2 and twist-3 two 
parton contributions for B — > ir form factors has been given in [491 ] . The factorizability of 
B — > 7T form factors under C* is shown valid upto twist-3 order under two parton approx- 
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imation. An all order proof of the factorization theorem for B — > tt form factors is shown 
valid at twist- 2 order in SCET [5^]. The factorizability of the twist-3 three parton contri- 
butions to the B — > 71 form factors is still inaccessible. We assume that this is also valid 
in this work. A simple analysis shows that the twist-3 three parton contributions are also 
factorizable in C* . However, the complete analysis is tedious and very technical. We skip 
this part and leave it for other place. In next section, the twist-3 three parton contributions 
from the pion are calculated based on this assumption. The analysis of the factorizability 
for the subleading twist contributions from the B meson is more complicate. The complete 
0(a s ) factorization analysis also needs a lot of space and is not appropriate to be given here. 
Similarly, we also assume that the factorization is also valid for the subleading twist con- 
tributions from the B meson. The calculations of the relevant quantities are based on this 
assumption. At order 0(a s ), the explicit results given in the next section for the considered 
contributions confirm the factorizability. However, a complete analysis for the subleading 
twist contributions is important to make sure that the calculations given in this paper are 
perturbatively meaningful. 

In Eq. (115p . the loop parton momenta of the partons of the B meson, Ib, is defined to 
flow from the B mesons into the scattering center, and those momenta I and I' are defined 
to flow from the scattering center into the tt meson. To perform collinear expansion, the 
momenta, I, Ib are found convenient to separate their on-shell part from their off-shell 
parts. For example, I can be written as 

I 2 + I 2 

p = P + l -±±Ln? + V\ (19) 
2n-l 

where I 2 denotes the virtuality of I, l^_ is the transverse momentum, and lj_lj_ )fJ , = —l\- 
I 2 = n 2 = 0. n is an auxiliary light-like vector. Under C* , only collinear partons can involve 
in hard scattering center. The most important contributions should come from collinear 



momentum configuration according to the power counting rules |41l . |43| . If I represents 
the momentum of a collinear parton, i.e., I is a collinear momentum, then / has a limited 
(transverse) virtuality, I 2 ~ l\ ~ 0(A 2 ). In C*, the partons external to the hard scattering 
center are on-shell. According to Eq. (fl9|) . the on-shell parton can have a momentum P with 
l 2 = 0, or, 

ll = P + + ll , (20) 

11 = 0. 
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The / M is the collinear component of / M and is used for collinear partons in . The lp has 
non-collinear momenta and is usually used in the kip. It is noted that n ■ I = uE with u the 
momentum fraction and E the energy scale. For u ~ 0(1), the n M term of lp is suppressed 
than the V 1 and V]_ terms. For u ~ 0(A/E), three parts of lp are equally important. There 
requires some cares for the use of lp in the kip. Especially, it is important in the end-point 
region. We argue that the n M term of lp should be kept in kip. The calculations for B — > -n 
form factors in k T should be rechecked. 

The CE is able to derive higher twist contributions from higher Fock state, non-collinear 
components of parton momentum, and miss-matched spin states. The higher twist con- 
tributions can be factorized into partonic and hadronic parts, which are separately gauge 
invariant. The parton model interpretation for higher twist contributions is similar to the 
leading twist contributions. 



The CE may confront with loop expansion (LE). We choose our strategy 4jJ, |43j to firstly 
perform the LE for the amplitudes and then to use CE. Changing the application order of 
these two expansions would not make differences. The parton functions are expanded in LE 
as 

H»(l Bl l) = H^(l B ,l) + 0(a 2 s ), (21) 
HZ(l B , I, I') = HW«(l B , I, I') + 0(a 2 s ) , (22) 

where and f/^ are of order 0(a s ). In this paper, we only consider 0(a s ) contribu- 

tions and assume 0(a 2 ) not important. 

The first step is to expand H^ ,tM (l B} l) and Hv'^(l B ,l,l') with respect to the collinear 
momenta, l B , I, and P, by a Taylor expansion 

arr(l),fi 

H^(l B ,l) = H^(U)+ -Q^\{ k =k, k =i B ,i}(l-iy + ---, (23) 

k=l,ls 

H^(l B ,l,l>) = H^(l B ,l,l')+J2 -^rl{^ fl ,M'}( fc -^ + --- > ( 24 ) 

k=l,V 

The choice of l B = n-l B is according to = n ■ to make the internal virtual gluon carry 
most violate energy. The Taylor expansion follows from the argument of the C* that the 
hard scattering functions contain only collinear momenta of partons (4]]]. The co- variant 
gauge, d ■ A = 0, is used. Up-to twist-3 order, one only needs to consider the first term of 
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the above first expansion series and the first two terms in the above second expansion series. 
By substituting them into the original convolution integrations, one can obtain 

M " = I duT ^ H(1) '"^ u ^)Mu)] (25) 
f dij duj du'Ti[H^{i,uy)MiH^ u {u,u'^ 



+ 
+ 



where ufl, = gZ — n v rn v . Two light-like auxiliary vectors are used = (n + ,n~ ,n±) = 

(1,0,0), n M = (n + ,n~,n±) = (0,1,0). The hard scattering functions H^^(^,u), 
H^'^^jUjU 1 ), Hl^'^(^,u,u') are defined by low energy theorems 

hW"(Z,u) = H^(l B ,i), (26) 

H^,u,u') = H^(i B ,l,i'), (27) 

H^u,v>) = _ a ^)\ {lB = Mv = P) . (28) 
The meson functions are defined as 

<MO = j^l d 4 z5U-^)e^(0\q(z)b(0)\B), (29) 

Mu) = Jj^J d*z5(u-^y l -(7r\q(z) q (0)\0), (30) 

(3l) 

x5(u - ^)5(u' - l ^)(n\q(z)(-g)A»(z')q(0)\0) , (32) 

Pw Pit 

x5(u - ^)5(u> - l -^)(n\q(z)(-ig)G^(z') q (0)\0) . (34) 

pi pi 



By using the identities 



,( tt -g) = /^*>, (35) 
f ^ f *** U ' { '-* n) = ( 36 ) 
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and their similarities for the corresponding integrations, the meson functions become 

MO = J ^e"€<0|?(An)6(0)|B>, (37) 

€M = /^^e tt '^<7r|g(^n)(-yK(^n) g (0)|0>, (39) 

<«>X) = f-^jf ?|^ eaVA ' (U '^ ) ^l?(^)(-^)^(^)?( )l )- ( 4 °) 

where = d'M 1 ' — d u A v . We note that the coupling g is absorbed into the meson functions, 
0^ and </>^. The three parton contributions are then counted as the same order of the two 
parton ones. This is different from the counting rule in 31], in which the three parton 
contributions are counted as one more 0(a s ) order than the two parton ones. 
The second term in Eq. (l25|) is related to gauge phase factors 

Jd^Jduj du'lk[H^{^u,u')Mi)€{^n')) (41) 

= / dif du'Ti[H^(^u')MO^,n-A^')}- J d£ J duTi[H^^,u)MO^,n-A^)} 
+0(A ± ) , 

where 

—e^(7r\q(-n)(-zg)j^ d V n ■ A a (^n)T a q(0)\0) . (42) 

0(A±) denote the terms composed of gauge fields with a transversal polarization and are 
identified as sub-leading twist contributions. The gauge phase factors are absorbed into the 
first term in Eq.( l25l) . In covariant gauge, there are infinite number of similar gauge phase 
factor terms from higher order Feynman diagrams. Their treatments are similar to the above 
and skipped here. The result becomes 

MM = 1^1 du ^ Hil) *te> u )M$Mu)] ( 43 ) 
+ jdijduj du'TT[Hl 1 ^^,u,u')MO^ u (u,u')} + --- . 

The color index factorization are performed in the following way, 

Tt c [H (1 ^^u)M04>M] = [H^^,u)} ijiM [MO}ik[M^hi, (44) 
Tr c [Hj»"(t, u, u')<f> B ^H<f>t v (u, u')] = u, u'^UMOhk^ [4>t >\u)] jm {^) 
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where i, j, k, I are color indices in the fundamental representation and b is the color index 
in the adjoint presentation. u)]ij t ki an d [H$'^' b (£, u, u')]^^ are expanded in terms 

of color factors 



(46) 



[H^>\t, u, u%j, km = H^(t u, v!) (J^5 ij 5 M (T b ) lm + (T%(T<%(T fc ) /m + ■ ■ ■ ) (.47) 

For 0(a s ) Feynman diagrams depicted in Fig. [H only the second terms in the right hand 
side of Eq. (I46|47D can contribute. They can be simplified by the color algebra 



1 Cp 
(T a )ij(T a ) k i(T b )i m = 5ik5 j i—Tr[T a T a ](T b )i rn = 5i k 5 j i — (T b )i m 



Eq. ((3SD is applied to Eq. (jUJ) to have 



SikSji—[4>B(0}ik[(P(u)}ji = —<Pb(04>M . 



Similarly Eq. (j4"9l is applied to Eq. fj45l) to have 

/rr b\ 



ikVjr 



c F 



(48) 
(49) 

(50) 

(51) 



Before the spin index factorization is performed, it is necessary to eliminate all terms 
in the H functions which may lead to higher twist contributions under the equation of 
motion for the quark. The equations of motion of light quarks (assumed mass-less), whose 
momentum is Z M in the collinear region Z M = ~ (E,A/E,A), are equivalently to 

the following identities 



P + 



ie 



(I - 1) ( l7a )_ —J— - ft 



l 2 + te 



I 2 + ie 

Since the propagator (the special propagator) 

i 



2n ■ I + ie 

i ri 



2n ■ I + ie 



i 1 



(52) 



2n • I + ie 

does not propagate, its associated terms are absorbed into the hard scattering functions (the 
last square bracket term in the last term of Eq. fl52|) ). The meanings of the above identities 
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are as follows. If there is one factor ft in the hard scattering functions contact with the 
long distance part of the quark parton propagator of the it meson, then the result is to 
extract one short distance part of the propagator and a vertex i^ a with a non-collinear 
momentum factor (/ — l) a . The non-collinear momentum factor will be absorbed by the n 
meson function, (j> n (u), to have 

(/ - irUu) = <> f^f ^^e^'-«)(7r|g(|A)^"'(|n)g(0)|0) (53) 



v. 



(54) 



The similar case arises when the ft in the hard scattering functions contact with the short 
distance part of the quark parton propagator of the n meson, the result is 



i ft 



2n ■ I + ie 



iy 



I' 2 + ie 



{-gA a ){i lc 



i t/i 



2n • I + ie 



(55) 



-(-gA a ') 



(Ha) 



2n ■ I + ie 



-n 



+ 



_l' 2 + ie 

where dots denotes the term would be absorbed by the gauge phase factor of the matrix 
element. The A a is the gauge field and its associated terms are defined to be absorbed into 
their corresponding hadron functions 

wZ(-gA°)M*) =u «f^f ^ elAvr?(M '" M) ^l^| A )(-^^(|^)^ )l°) 

= <,C>X) . 

The total effect, when one fl factor can contact with the partons of the 7r meson, is 



Tr [#«(£, tt)<MO0J 



u 



(56) 



Tr 



2n ■ l q + ie 



2n ■ lg — ie 



-< 7 «) ^Mm^H 



where 



<o(u,u>) = J^I g^e^'-«)(vr|g(|A)^'(|n)g(0)|0), (57) 

where iD a ' = id a ' — gA a ' . Since n ■ l q ^ are of order E for collinear l q ^, the related 
contributions are suppressed by one additional i? _1 order. We call these terms as abnormal 
terms. The other terms are identified as normal terms. The normal terms will be kept in 
the reduced hard scattering functions, while the abnormal terms are dropped. 
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In the form factors, the £ is of order 0(A/E) due to the twist- 2 B meson distribution 
amplitude, 0^" 2 (O- The expansion into short distance and long parts of the quark propagator 
can not help to separate different twist contributions from the B meson. In this work, we 
only consider the contributions from the two parton Fock state \bq) of the B meson. The 
contributions from the sub-leading twist state of B meson are left to other places. The 
equation of motion for the b quark, (ft, — m^b = 0, is the only condition. This fact reflects 
in the derivation of the reduced hard scattering functions. (Refer to following text and 
Appendix A.) 

Considering all possibilities of the applications of equations of motion of quarks, the 
spin structures of the hard scattering functions are strongly restricted. It is useful to take 
the diagram in Fig. 1(a) as an example to explain this operation. After the color index 
factorization, the amplitude for Fig. 1(a) is proportional to 

J dij duTr[HW^ la \Z,u)cj> B (0<t>M} ■ (58) 
The spin structure of H^'^ 1 ^ (£,u) has the expression 

[H^ la \C,u)] ijM oc WCK- #h%ha]ki, (59) 
where ]/ q = u fa and = ifin ■ P h — u fa. By using 

g ap = n<* n P + n "nP + df , (60) 
7 « = n a i + n a fL + -il , (61) 

where 7" = (^73 and g% — d™ — —2. The spin structure of if( 1 )'^'( la )(£ ) u) is then expanded 
into a series in terms of ft, /i, 7°. Each term in the expansion series is then examined to 
determine whether it is of the considered twist or of higher twist according to the equations of 
motion for the quarks in the n meson and the B meson. The abnormal terms are subtracted 
from the expression. We note that this analysis is automatically fulfilled for leading twist 
hard scattering functions, because the leading twist spin structure of the 7r meson can 
eliminate the possible fi factor by fa ft ocfi ft — 0. For the sub-leading twist amplitudes, 
this procedure of subtraction of possible abnormal terms is necessary. For example, the spin 
structures of the twist-3 PS or PT DA are proportional to 75 or e^^a^n^vP . These spin 
factors can not eliminate the abnormal terms by the traditional method. In the following, 
we assume that the hard scattering functions are determined according to the above analysis 
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and the resultant hard scattering functions are called reduced hard scattering functions. The 
relevant reduced hard scattering functions are given in Appendix A. 
The spin index factorization is performed in the following way, 

Tv s [H^(^u)MOMn)} = [H^(Z,u)}^ M [MOWMu)h, (62) 
Tr s [H^,u,u')MO^ V (^^ = W ) ^(e,«,« , )]«,H^B(0]*<[^ , '(«)]ii, (63) 
where the meson functions are expanded as 

[MZ)]ik = ^(7 P 75kTr[<M07p7 5 ] + ^TsMMOts] + " ' " > ( 64 ) 
[M»)] fl = \{l P l,)ji^[4>MlM + \{l^[<t>,{0^] (65) 

+\{.Vp\ib)jf£A ( t ) M(j px ^} + ■■■ , 

'<t>t{u,u') 

The dots denote those terms are not of our interesting. Each coefficient in the above expan- 



sion is attributed by a DA according to the following definitions 

Tr[<M07 P 7 5 ] = if B [P B , P MO + E{n p - n p )^ B ] + • • • , (67) 

Tr[<M07 5 ] = if B m B MO + " • • , (68) 

Ti-[^ W («)7 P 7 5 ] = -iUPn, P ^{u) + • • • , (69) 

Tr[0.(O7 5 ] = + (70) 

TrlMu)*^] = -f n » x [fi p ,n x ]<%(u) + ■■■, (71) 

Tr[<f>t u (u, u')a px l5 ] = -/^(tfrfV - p x d^)p^(u, «') + •••■ (72) 



In these definitions for the DAs, only relevant ones are shown and the others are left into the 
dots terms. Note that the 0^ ,ct (m) are defined as the the energetic limits of the PS and PT 
DAs of the n meson, respectively. The spin projector of is the leading part of the full 

spin projector [p£, z x ] of the PT DA under the energetic limit. The (/>°(u) corresponds to the 
d(f)°' c (u)/du of the usually used PT DA, <J)^ c (u). The factor /j, x = m^j{m q + m g -) with 
the pion mass and m q ^ the current quark(anti-quark) masses. The Dirac matrices in the 
expansion series are absorbed by the hard scattering functions. The spin index factorization 
is completed. 

The final result is written as 

J^fJ. = jy[tw2,n _|_ j^B,p _|_ j^tw3,ps,p _|_ j^tw3,pt,fi _|_ ]^-tw3,3p,n _)_... (73) 
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il 



-(a pXl5 ) jl Tr[<f>^(u,u')a? x l5 ] + 



(66) 



where the expansion terms in the right hand side of Eq. (!73|) are defined as 

M tw2 ^ = ^fsUj d£MO J duj>$(u)Tr [75 ^H^(^u)(^ B + m B ) l5 ] , (74) 

M B >» = -^f B UE J d£MO I duti(u)Tr [75 ^H^{^u){fi- fi) lh \ , (75) 

M tw3 > ps ^ = ^fBUHxf d£MZ) I durM^hH^(^ U )(F B + m B ) l5 ] , (76) 

M t W 3, P t,» = _±_ fBf ^ x j d ^ B ^ J durMTr[e ± -aH^(Z,u)(F B + m B ) l5 ] , (77) 

M W P ,» = _l_ fBf ^ x j d ^ B ^ J dUq j du _ j dUg ^p{u q)U ^ Ug) (78) 
X8(l -u q - u- q - i£ 9 )Tr [ W5 #«^(£, u q , u q , u g )(?> B + m B ) 75 ] wjF^foO , 

where r" A " v (p 7r ) = (p^d^' -pld^')p v K and e± ■ a = e ± ^a al3 = e^a^vPv? with a a ? = 
i[ 7 a )7 ^]/2. M tw2 ^ has been calculated in previous section. The other four twist-3 terms, 
jyj-tw3,B,ps,pt,3p,n^ are ca i cu i a ted in next Section. In the above expression, we have used a 
transformation of the integrations over u and u' into the integrations over u q , u q , and u g to 
make the expression more explicitly. The u q , u q , and u g are defined as the fractions of the 
quark, anti-quark, and gluon partons of the 7r meson. 

The retain of £ in the internal b quark propagator can partially resolve the end-point 
divergences at twist-3, which are linear divergences due to the constant nature of the ir 
meson's twist-3 distribution amplitude (the pseudo-scalar one, 0^(w)). The PS DA 0^(w) = 1 
for the 7r meson is determined by the equation of motion for the quark of the 7r meson in 
the soft pion state, where the energy of the pion, E n , the pion mass, m n , or the scale of the 
parton transverse momentum, Zj_, are of the same order, E n ~ m n ~ l± ~ 0(A) . The final 
state 7r meson in the semi-leptonic B° — > tt~1 + vi process can carry an energetic momentum, 
where E n 3> m^. The equation of motion for <fiP(u) can be further approximated to a reduced 
equation of motion 4l|. According to this reduced equation of motion, the it meson's PS 
DA and the PT DA, <fi°(u), become equal. In this way, 4>^(u) and 4>°{u) can be modeled 
similar to the twist-2 tt meson's distribution amplitude, <p^(u), to take their asymptotic 
forms, <fi^{u) = (p^iu) = 6u(l — u). The constant feature of 0^(w) as found in literature is 
disappeared in the n meson's energetic state. 

In the 7r meson's soft pion state, different twist order contributions ordered by 0((A/i? 7r ) n ) 
become equally important under the limit E n ~ A. On the other hand, in the 7r meson's 
energetic state [41], different twist order contributions ordered by 0((A/E n ) n ) are restric- 
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tively ordered under the limit E n ^> A. This explains why </j^(u) can have different forms 
in the energetic and chiral limits. Once the non-constant <f)P(u) are substituted into the 
relevant expressions for the amplitude under the (, R , the linear divergences automatically 
disappear and the related end-point divergent problem is resolved. There are no similar 
end-point divergences for the twist-3 three parton amplitude M tw3,3p ' 11 since the three par- 
ton DA (p^iuq.Uq^Ug) is not a constant. As a result, we obtain a factorization formula for 
the B — > 7r form factors up-to twist-3 order. 
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Figure 3: Feynman diagrams for the hard scattering functions of the C* amplitudes for the 
B — > ttW[ decays. The external meson states \B) and (ir\ are not shown. The cross vertex denotes 
the vector operator 7^. The dot vertex denotes the subscript 77 in H^„. Only contributions from 
diagrams (1-4) are of twist-3 order. The contributions from diagram (5-36) are of higher than 
twist- 3 order. 



The calculations of three parton contributions are straightforward by using the reduced 
hard scattering functions given in Appendix A. The important steps in the calculations 
are described below. There are totally 36 Feynman diagrams as depicted in Fig. [3j Only 
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4 Feynman diagrams (Fig. 0(1)- (4)) are needed at twist-3 order. The rest 32 Feynman 
diagrams (Fig. E(5)-(36)) are of higher twist. Under covariant gauge, the contraction of the 
spin projector T pXriu (p n ) = (p^d^ — p x d pv )p^ with those contributions from Fig. E](l)-(4) 



contained in Tr 



J tig, Uq^ llg 

)(f> B + m B )l5\ results in Eq.(l9Tj) 



IV. TWIST-3 AMPLITUDES AND FORM FACTORS 

In this Section, we present the twist-3 amplitudes M B ^, M tw3 ^, M tw3 > pt ^, and M tw3 ' 3p > p 
calculated by the collinear expansion introduced in last Section. The twist-3 form factors 
/*™ 3 are extracted from each twist-3 amplitude. The twist-3 form factors p Bl ^ tw3 are d e _ 
termined from the twist-3 form factors fffi. The calculations of Mf? 3 * 3 , M B '^, M tw3 ' ps '^, 
M tw3 ' pt,p , and M tw3 ' 3p,p are straightforward according to Eqs. (17611771179]) . The reduced two 
parton hard scattering functions H^' p (^,u) and the reduced three parton hard scattering 
functions H^'^ are given in Appendix A. By substituting these reduced hard scattering 
functions into Eqs. (l76ll77fT9l) . one can obtain the following twist-3 contributions. 



A. (pB contributions 



The 4>b contributions are written as 



M 



Mb fdfMO f du^(u)H B ^, 
Jo Jo 



where 



e b^ = 1 p n , (l + y)yu-T]- (1 -r])j ^ 



Similar to the twist- 2 case, the form factors are defined as 



(79) 



(80) 



where f B 2 are given by 

ira s Cp 

ira s Cp 



(0 J du<p P {u) 



(1 + T])7]U - 7] - (1 -T})£' 

r) 2 £u(£ — t)u) 



(81) 
(82) 
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The form factors F+ Q (q 2 ) are obtained by f{ w 2 3 ' B as 

7TCt s Cf 



W) = 



2E 2 N, 



rj u-rj-(l- 2r])f 
r] 2 £u(£ - rju) 

£- (1 - 2rj)7]u - rf 
ri£u(Z - rju) 



(83) 



(84) 



B. Twist-3 two parton pseudo-scalar and pseudo-tensor contributions 

The pseudo-scalar contributions and pseudo-tensor contributions are found identical as 

'KOL s Cf 



j^tw3,ps,^i _ j^tw3,pt,ti = f f B 

ZJ\ C 

where 

The form factors are defined as 



27V, r * 



1 



Jo Jo 



jy-tw3,ps(pt),n _ ftv>3,ps(pt)^2jpg _|_ jtw3,ps(pt) 

where the form factors y*™ 3 > ps ( p ') are given by (2E 2 = m 2 B and r x = 2ji x /m B ) 



n 



tw3,ps(pt) / 2 



f: 



tw3,ps(j>t) 



na s Cp 
'2E 2 iVc 

7ra s Cf 
'2^2 W 



Uf B r x j d£MZ) f du<j> P M 



The form factors F+Q ,tw3,ps (q 2 ) are obtained by f^" 2 ^ VS as 



7 B-K,tw3,ps(pt) (n 2 



(q 2 



na s Cp 



i(0 / 



u 



(1 +r])r]u-£ 
r] 3 £u(£ - rju) 



rju- (1 — r/)€ 
rf £w(£ - rju) 



B-K,tw3,ps[pt), 2\ M s C F 



V) 



4E N t 

1 - 2r])r]u - (1 - /?)£ 



x 



(85) 



(86) 
(87) 



(88) 



Uf B r x j d^ B (0 I dwft(u) (89) 



ri 2 £u(£ - rju) 

We found that the twist-3 two pseudo-tensor amplitude M twZ > pt ^ is equal to the twist-3 
two pseudo-scalar amplitude M tw3 ' ps '^. The corresponding form factors defined by these two 
amplitudes are identical. This is consistent with the reduced equations of motion (fiP(u) = 

#(«)■ 



25 



C. Twist-3 three parton contributions 



The twist-3 three parton contributions are written as 

ttq: s Cf 



N r . 



-A/b^ / d£<f) B {£) j du q J dug J du g 



X5(l-U q -Uq- U g )(j)l P (u q , Uq, U g )H 3p U q , Uq, U g ) , 



where 



H^,Uq,Uq,U g ) 



1 



tlqllqllg 



^rfugE 2 



(1 + r]u g )r]Ug -£ 



PS 



The twist-3 three parton form factors are defined as 



where the form factors /f^ 3 ' p are given by (2E 2 = m 2 B and r x = 2fi x /m B ) 



/f 3 ' 3 V) 



7TQ; S Cp 



E 2 N t 
4>l P {u q ,Uq, Ug) 



tlqllqllg 



7] 2 ^Ug 



7ra s Cf 



fnfBr x J d£(p B (£) J duq j dug j du g 8(l 



f(Ug,Uq, U g ) 



(1 + r]Ug)r]Ug - ( 
V 3 &g(£-r]Uq) 

The form factors F^ tw3,3p [q 2 ) are obtained by Eq.(j3 



UqUqUg 



-^Bir,tw3,3p ( n 2 
+ 



(O 



2E 2 W r 



(90) 



(91) 



fnfBr x J d£<j> B (£) J du q J duq j du g 5(l - u q - u q - u g ) (92) 



Uq-Uq- Ug) (93) 



fwfBrl J d£4> B (£) J du q j duq j du g 5(l - u q - u q - u g ) (94) 
(1 - rpig)TfU q - fj£" 



7T \ u q, Uq, Ug 
UqUqUg 



V 3 £u g {£ - 7]Uq) 



-iB-n,tw3,3p ( 2 



nct s Cf 

Ie 2 ^. 



USbT^ j d£<AB(£) j dUq j dUq j du g 5(l - Uq - Uq - Ug) (95) 

fj£- (1-7/(2 -u g ))r]u q 



71- V^9' ^g> 
UqUqUg 



V 2 £Ug{£ ~ TJUg) 
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D. The resultant form factors 



The form factors up-to 0(1/ E) are 



F»*(q 2 ) = F t + w2 (q 2 ) + F»(q 2 ) + F^ ps (q 2 ) + F™ 6 ' pt (q 2 ) + F t + W ^ p (q 2 ) 



-itu)3,pi/ 2\ 



(96) 



na s (t) Cp 



x 



2E 2 N t 

ixa s (t) Cp 
2E 2 W c 

7ra s (t) Cp 
2E 2 W c 

(f>l P (u q ,U q , U 



/tt/b 



^( u ) + + <K(u)r^[(l - r])£ - rju] 



/tt/b / d&B 



(f) / 



ducj)^(u) 



r] 2 ^u(i - T)u) 

rju — i] — (1 — 2rj)^ 
r] 2 £u(£ - rju) 



UfBrZ J df<MO J du q J du q J du g 5(l 
(1 - r]u g )vuq - vC 



U q - Uq 



UqUqllg 



V 3 £u g (£ - r)Uq){, 



F^(q 2 ) = F* w2 (q 2 ) + F B (q 2 ) + F^' ps (q 2 ) + F^' pt (q 2 ) + F to3 ' 3p (g 2 ) 
na s (t) Cp 



(97) 



/tt/b / d£(f) B 



X 



2E 2 N c 

[77 + (1 - 2tj)£ - (1 - 2^)^] + 0P( M )r x [(1 - 277)77^ - (1 - V )£\ 



ira s (t) Cp 
2E 2 iVT 
7ra s (t) Cf 

2£ 2 iv: 



rfiud - r)u) 

£ — (1 — 2r))r)u — i] 



fnfBr* I d£(p B (0 I duq I duq 



f]iu(i - rju) 



du g 5(l — U q — Uq — U g ) 



X 



7{U q ,Ug, U g ) 

UqUqllg 



v£ - ( x - v( 2 - u g ))yu q 
V 2 £u g (£ - r)Ug) 



The form factor F^' K (q 2 ) is analytic in the whole range of q 2 , < q 2 < q max , g^ax = 26.42 
GeV 2 for m B = 5.28 GeV and m w = 0.14 GeV. Under the maximal recoil limit, rj max = 1, 
the form factor F? n becomes 



Inn^fV) = ^%Ub J dZMO I du 



£(«)£+ #(u)(l-rj)ii 



na s (t) Cp 
2E 2 iVT 
7ra s (t) 

2£ 2 ivT 



(98) 



Ub J d^MO J du^(u) 

/tt/s^ j d£(f>B(0 J du q J dUg J du g 5(l - U q - Ug - U g ) 



X 



7T V M <2> M <2> u ffJ 



UqUg 



1 



The dependence of strong coupling constant ai s (i) on i oc ^fr\ restrict the effective range of 
q 2 to be < q 2 < 16 GeV 2 , in which the values of a s are varying within 0.33 ~ 0.48. 
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V. NUMERICAL CALCULATIONS 



We apply previous results for numerical calculations of the form factor F+ n (q 2 ) for B° — > 
7r~/ + z// decays. The shape parameter ujb of the B meson DA and the value of \V u b\ will be 
determined. The q 2 shape of the form factor F^ n (q 2 ) is given according to the value of ujb- 



A. analysis method 

In order to explore the full range q 2 shape of F + (q 2 ) based on our calculations, we employ 
a statistic analysis method. We assume that the form factor F + (q 2 ) can smoothly vary 
with q 2 from small q 2 region to large q 2 region. We first collect 17 data points of F + (q 2 ) at 
< q 2 < 16 GeV 2 and Aq 2 = lGeV 2 . We denote these 17 data points as the C f set of data 
points. Based on the set of data points, we derive a fit formula for F + (q 2 ) by a minimum 
X 2 fit according to a Xff function. The fit form factor F^(q 2 ) is a function of q 2 with 9 
parameters, Oj, i — 0, 1, - • • ,8, as 

8 

Ff(x) = ^2a i x i . (99) 



i=0 



The Xff function is defined as 



j^ (Fl\q 2 ) - F^(q 2 )) 2 ^ (10Q) 



Xff , , 

o 



where % denotes the i-th bin and N is the total bin number of data points, q 2 and Oi are 
the reference q 2 value and the uncertainty for the i-th bin of the used data set. F^(qf) is 
the value of F+ lt (x) at x — q 2 . F^ ata (g 2 ) is the z-bin data. A 8% error is taken for the input 
parameters (a s , f n , Ib, r x ). See more detailed explanations about this 8% error in next 
Subsection. 

The correctness of the q 2 shape of F^(q 2 ) for q 2 > 16 GeV 2 is estimated by comparing the 
oredictions of -Fj^g 2 ) for q 2 > 16 GeV 2 to the LQCD calculations by the FNAL collaboration 
Q] and by the HPQCD collaboration [l^], respectively. 



16 



The fit form factors is used to make a prediction for the branching ratio denoted as 
Br§l for B — ► -nlv\. The ujb and \V u b\ are determined by comparing the predictions and 
experimental data, Br^. This is equivalent to a fit for the value of |V^b-P+(0)|. Because 
the q 2 shape of F + (q 2 ) is fixed by Eq. (1961) . theoretically, the normalization point |V^&F + (0)| 
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can be determined by the fitting method. The |F + (0)| and \ V u b\ are determined separately. 
Our best fit to the central value of Br e ^' ful1 = (1.36 ± 0.09) x 10~ 4 (PDG08) [45J gives 



\V ub F+{0)\ = (6.471 ± 0.196 cxp ± 0.523 th ) x 10" 4 , 

|F+(0)| = 0.164lS| a .±0.008| /fl (101) 
= 0.164±0.013, 
\V ub \ = (3.946^;H cxp ^l th ) x 10- 3 , 
lo b = 0.315GeV. 

The errors are experimental (exp) and theoretical (th). The theoretical errors are only 
considering those uncertainties from a s and |F+(0)| is calculated according to the fit 
form factor and the fitted value ujb- \V ub \ is determined by a combination of |\4b-^+(0)| and 
|F + (0)|. The central value of |146-F + (0)| is determined by a best fit to the central value of the 
Br c ^ . The determinations of the fit form factor and corresponding \ V u b\ will be respectively 
described in detail latter. 

The form factor F^ f (q 2 ) is perturbatively meaningful for q 2 < 16 GeV 2 . We perform 
an similar analysis as above by using the partial branching fraction Br^' q <16GcV = (0.93 ± 
0.06) x 10~ 4 . 1 The best fit result is given by 

\V ub F + {Q>)\ = (7.34±0.233 cxp ±0.597 th ) x 10" 4 , (102) 
|F+(0)| = 0.188^oU s ±0.009| /s 



0.188±0.014 



th 



\V ub \ = (3.904±g;iH c xp ^ th )xl0- 3 , 
co B = 0.281GeV, 

The errors are the same as explained previous. We note that our best fit value \V ub F+(0) \ 



6.5 x 10 4 is lower than the value 23J], which was obtained by a fitting to the branchin 



fraction and q 2 spectrum data by BaBar collaboration 53 ] 



VriFfT (0)| = (9.1 ±0.3 ±0.6) x 10" 



1 The value is quoted from the ICHEP08 averages given in the online update at 
|http : / / www .slac.stanford.edu/xo rg/hfag by the HFAG colaboration. Note that the full branching 
fraction Br e s ^ Ml = (1.34 ± 0.08) x 10" 4 has been updated by HFAG at 2008 
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where the first error is from the uncertainty on BriB — > ttIui), and the second from the 
parametrization of shape of the form factor versus q 2 . The |^bF + (0)| in Eqs. fllOip and 
PUD is consistent with the result |^F + (0)| = (7.6 ± 1.9) x KT 4 obtained with SCET and 
QCD factorization j^]. 

The analysis procedure is the following. We first choose a reference value of uj b to prepare 
the G* set of data point for F + . Using this C* set of data, the fit form factor is determined 
according to the Xff function. The fit form factor is then used to make a prediction for 
the branching fraction for B — > ixli>i denoted as Br<£ by varying the value of \V u b\. A best 
fit for the predicted branching fraction and experimental one is performed according to the 
XB r function 



( Br fit - Rr exp ) 2 

Xlr = ( SL 2 (103) 



where a is the experimental error calculated by adding the systematical and statistical errors 
in quadrature. Once a best fit is found, \V u b\ is determined. Otherwise, another u B is chosen 
and the whole procedure is run from start, again. 



B. form factors 



The values of a s (t) at different values of q 2 are calculated by the program [54J and given 
in Table [H The program 541] has considered experimental and theoretical uncertainties. The 
scale t is defined as the scale of the internal virtual gluon, which are assumed to carry a 
hard-collinear energy, t = \J r i]h} l m B with = 0.5 GeV and m# = 5.28 GeV. 

The uncertainties on the F + could be (1) phenomenological models for meson DAs, (2) 
input parameters (a s , f n , r£, u B , r] 3 , and u>). The uncertainty on phenomenological 
models is fixed by taking models determined from other processes. The twist-2 and twist-3 
meson DAs takes the following models 



#(«) = ^(«)=^(«) = 6u(l-«), 



UJ 



w p (u q ,Uq,u g ) = 360i] 3 u q UgU g (l + -(7u g 



3)) 



where rj3 = 0.015 and u = —3.0 are used. The chiral factor is set r£ = 0.76 GeV. The 
total uncertainties on the input parameters are assumed constant of 8%. The recent LQCD 
calculations for f B are f B = MeV (CP-PACS) [55], f B — 190 i 14 MeV (JL-QCD) 
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Table I: The values of a s at q 2 , 0< q 2 < 16GeV 2 , are calculated by the program [54| 



g 2 (GeV 2 ) 





2 


4 


6 


8 


10 


12 


14 


16 


a s (t) 


0.333 


0.340 


0.348 


0.357 


0.368 


0.400 


0.421 


0.447 


0.484 


6a s (t) 


+0.020 
-0.018 


+0.20 
-0.019 


+0.022 
-0.019 


+0.023 
-0.021 


+0.025 
-0.022 


+0.032 
-0.028 


+0.036 
-0.031 


+0.043 
-0.036 


+0.053 
-0.043 



54 f B = 206 ±10 MeV (ALPHA) [57], f B = 216 ±11 MeV (HPQCD) |58|,|59|, f B = 214 ±9 
MeV (Guo-Weng) 6(3]. Except of the CP-PACS result with 16% errors, the other calculations 
contain about 5% uncertainties. In our calculations, we choose f B = 200 MeV as a reference 
value and associate a 5% error on it. The chiral enhanced contributions are about 20% in 
our calculations. We choose r£ = 2m 2 K / {m B {m u ±m<f)) = 0.76 GeV. This corresponds to use 



451 ] . rriir = 140 MeV, and m B = 5.28 GeV. The pion decay constant 



(m u + m d ) = 9.8 MeV 

,/U = (130.4 ± 0.04 ± 0.02) MeV and f n o = (130 ± 5) MeV 45]. The errors on f n - can be 



neglected. The uncertainties on a s are 5 — 10%. There are about 10% errors associated with 
r x . Because the chiral contributions with r x are about 20% of the leading twist (twist-2) 
contributions, the errors on r x are about 0.2 x 0.1 ~ 2% and can be neglected safely. The 
constant 8% errors on F + (q 2 ) is assumed according to the errors on a s and f B . Instead of 
using A = m B — rrib and X B for the value of uj b , we take uj b as a free parameter and propose 
to determine uj b and \V u b\ by the world branching fraction BrT? for B — > -nlv\ 



Br e s l P (B° -> u-tut) = (1.36 ± 0.06 stat ± 0.07 sys ) x 10" 



(104) 



for < q 2 < q 2 max , where q 2 max = 26.42 GeV 2 . 

17 bins of G* data points of F+ n (q 2 ) with Aq 2 = lGeV 2 are prepared by taking a 
test value of uo B . In our present case, we choose uo B = 0.46 GeV as the first value. The 
i+ (g 2 ) is then used to make a x 2 fit to these 17 bins of data. The best fit is determined 
by requiring the Xff function to be minimal. The fit form factor is then substituted in 



Eq. (I109j) to make a prediction for the branching fraction, Br$l by varying the value of \V V 



ab 



in 2.8 x 10 3 < \V u b\ < 4.5 x 10 3 . The best fit is found by using the xlr function to calculate 
the minimal \ 2 value. Our best fit gives 

uj b = 0.315 ± 0.000 GeV , \V ub \ 



3.946l°;™2 x 10- 



(105) 



under x| r < 0.001. The errors on the fit value of uo B and \V u b\ are statistic. The fit form 
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factor is assumed to have 9 parameters 



.-■fit, 



x) 



X 



(106) 



The fitting of Ff t (x) to the 17 data points is performed by using different numbers of 
parameters, 1 to 9. The minimal \ 2 value of the Xff function is required to derive a best fit. 
The Xff/N is equal to 0.05/ (17 — 5) by assuming an 8% constant errors on the C* data set. 
The values of parameters <2j, i = • • • 8, are given in Table [Til The errors on ai, i = • • • 8, 
are statistic due to the best fit of uo B . We denote the fit form factor as F^ \q 2 ). 

By F + -'(q 2 ), we compare our calculations with the same form factor calculated by the 
other approaches, including the kj, method 31], the LCSR method 



method 



13| the LQCD-Fermilab/MILC/2005 (FNAL05) method 



61 



16 



the LQCD-HPQCD 



5l| and the LQCD- 



Fermilab/MILC/2008 (FNAL08) method 52J. The comparisons are shown in Fig. HI The 
LCSR form factor 6l| is expressed as 



-tLCSR/^2 



III T 



+ 



r 2 



(107) 



where r x = 0.744, r 2 = -0.486, m l = 5.32 GeV, m| t = 40.73 GeV 2 . Both LCSR and C f 
form factors are larger than the LQCD form factors at most range of q 2 . The form factor 
is larger than the C* and LCSR form factors in the available range, q 2 < 10 GeV 2 . We find 
that the \ 2 /N values of Ff(q 2 ) versus the FNAL(=FNAL05+FNAL08) and HPQCD data 
are 54.8/15 and 271.1/7, respectively. It is seen that Ff(q 2 ) for q 2 > 16 GeV 2 is close to 
the FNAL data but has a large deviation with the HPQCD data. 

The twist-2, j?^ n ^ tw2 ^ anc i twist-3, p^ n ^ tw3 ^ components of the form factor F^ w are calcu- 
lated at different g 2 , < q 2 < 16 GeV 2 in Table. [TTT1 It is seen that F_ 



Bir,tw2 



> F_ 



Bir,tw3 1 



FBir,tw3,ps . r^Bn,tw3,pt . j-^Bir,tw3,3p , j-^Bir,tw4,B\ • ,1 -j j c 1 ti 

+ + r + + r + + r + ) in the considered range of q . lhe power 

expansion is meaningful. We note that F^ f (0) = 0.164 is obtained. This is close to the 

finding by the SCET analysis for B -> tttt decays [36], that F| CET (0) = 0.170. 

We note that F% f (0) = 0.164 is smaller than F+(0) = 0.25 ~ 0.30 derived by other 

theories (such as, LCSR, LQCD, fcj,,), and experiments. As was firstly found by Bauer et 



a/36j, there is an inconsistency between the small F + ET ' (0) = 0.17 implied in B — » nn 
decays and the large F + (0) = 0.25 ~ 0.30 derived by other theories and experiment. For 
example, F?*(0) = 0.28 was used in the QCD factorization for B — > nn decays and the 
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Table II: The values of parameters of F+ tt (q 2 ) and the best fit value of u>b in unit of GeV. The 



best fit means Xbv — 0-001. The errors on and a,i, i = 0, • • • ,8, are statistic. 



method 


2 

^ff- ujb x 10 1 a x 10 1 ai x 10 3 a 2 x 10 3 a 3 x 10 4 a 4 x 10 5 a 5 a 6 a 7 a 8 


full 


2j2p 3.151±| 1.636±| 8.403i}g 1.903±| -1.873±| 1.407±| 


partial q 2 < 16 GeV 2 


2f 2.810+J 1.881±| 9.773±|| 2.185±| -2.143±| 1.625±| 



LQCD(HPQCD) 
LgCDfFNAL05 
LQCD FNAL08) 

LCSP^ 



10 



15 



20 



25 



Figure 4: Comparisons between (solid line), kL (dash-dot-dash line), LCSR(dot line), , LQCD- 
FNAL (data point), and LQCD-HPQCD (data point) form factors . 



QCD factorization predictions for the branching ratios of B — > tttt decays were inconsis- 
tent with the experimental data. On the other hand, if F^ CET (0) = 0.17, then the QCD 
factorization based on SCET can explain the experimental data consistently. According 
to our calculations, the small F+ f (0) = 0.164 ± 0.013 can also explain the semi-leptonic 
decay B — > wWi with a value of \V u h\ consistent with the world averaged value. Recently, a 



general parameterization approach [14 



15 



16 



17 



has been used to st udy the q 2 shape 



of F 4 
et al 



q 2 ). Among of different parameterization approaches Q, 16, 17, 18], Arneson 
16] found that it is possible to consistently explain the branching ratio and BaBar 
differential rate data by using the SCET value F+ CET (0) = 0.17. The other parameteriza- 



tion approaches 



15 



17 



181 ] employed the LCSR value F + (0) = 0.26 as an input. Whether 
the value of F + (0) is small or large is a controversial topic, it needs more theoretical and 
experimental works to clarify this problem. 



33 



Table III: The values of different contributions of \F^(q 2 )\ at q 2 , 0< q 2 < 16GeV 2 . lo b 
0.315 GeV is used. 



g 2 (GeV 2 ) 


2 4 6 8 10 12 14 16 


\F^' tw2 (q 2 )\ x 10 3 

\F^ B (q 2 )\ x 10 4 
\ F Bir,tw3,ps+pt( q 2^ x 1Q 4 

|Ff r ' tu,3 '*( g 2 )| X 10 4 
total xlO 3 


175 205 240 286 351 444 582 787 1100 
129 152 180 218 270 347 463 640 919 
196 248 318 417 568 810 1213 1911 3202 
106 127 154 193 255 360 547 900 1615 
162 188 219 259 314 393 508 679 943 



The fit result u B = 0.315GeV implies A = 0.471 GeV or X B = 0.315GeV. We note 



that A s = 0.315GeV is lower than 0.460 ±0.110 GeV by Braun et a/|62j, 0.454 GeV by 



Grozin and Neubert 



46j, 0.479 ± 0.089 GeV by Lee and Neubert [47[, 0.600 GeV by Ball 
63j . nib = m B — A = (5.28 — 0.471) GeV = 4.57 GeV, is consistent with the m b determined 
by different schemes, m\ s ,= 4.701 ± 0.030 GeV (IS scheme)[64|, mf F = 4.630 ± 0.060 GeV 

4.613 ± 0.035 GeV (kinetic scheme)!^, but is 



641 ]. and 



m 



kin 



(shape function scheme) 

larger than m^ IS = 4.164 ± 0.025 GeV (MS bar scheme) 65jj. At the chiral point, E w ~ m- 



and q 2 (E n ) = 26.42 GeV 2 , the chiral perturbation theory predicts that 

9fsm B 



^Chiral (2 



(Q 2 



(108) 



2f 7T (E 7r + m B * -m B ) 

where g is the B*Bn coupling. The fit form factor F^(q 2 ) predicts F + * (26.42) = 5.12, 
which implies g = 0.234 ± 0.019 with a theoretical error by using m B * — m B = 45.75 ± 0.35 
MeV. The extracted q = 0.234 ± 0.019 is consistent with g = 0.22 ± 0.07 determined by the 



FNAL collaboration 52j and < g < 0.45 by the HPQCD collaboration 13], but is lower 
than the usually employed value g = 0.51 proposed by Stewart et al 



C. determination of V u b 

The differential decay rate for B° 
masses are vanishing, is given by 

dT 



7i under the approximation that the lepton 



dq 2 



th 



r 2 it/ i 2 

(j F\ v ub\ | ,3, wB -k, 2\\2 



(109) 
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where p n is the momentum of pion in the B meson rest frame. The branching ratio Br(B° 
it~l + ui) is expressed in terms of (dT/di])^ 



("Lax / dV\ 

Br(B° -> 7T-/+Z7,) = TB o J dq 2 ) 



) th 



(110) 



where r^o = (1.530 ± 0.009)ps is the life time of B° meson 45] 
To determine \V u b\, we employ 

/ dT \ 



R(q 2 v 



PQmax 



th 



The result is 



(5.709 ± O.giSth)^" 1 



By substituting R(qmax) Eq. (IllOp . we can determine 



\v: 



ub 



(3.946 ± 0.117 cxp ± 0.312th) x 10" 



111) 



(112) 



;ii3) 



The first error on \V u b\ is from branching ratio and the second error is from F + (0). The 
theoretical uncertainty is 8% in the same level with that of the inclusive method. For 
comparison, the theoretical errors are 15% in the LCSR method, 10 — 14% in the LQCD 



451 ]. and 12% in the parameterization approach (PA) 

ub 



14 



15 



(11131) is consistent with the inclusive value |V^ ntl 



avg 
ub 



B — > X u lui [45!] and the world averaged \V l 
the world averaged exclusive value is |V^ cl | 



The \Vub\ in Eq. 
(4.12 ± 0.15 cxp ± 0.40th) x 10" 3 from 
= (3.95 ± 0.35) x 10" 3 [45J. We note that 
(3.5 ± 0.6) x 10~ 3 45]. The discrepancy 



between the values of \V u i,\ determined from exclusive and inclusive methods has raised 
a lot of discussions in literature [23]. However, we note that a smaller inclusive value 
|VX 108 | = (3.98 ±0.15 ±0.30) x 10~ 3 with experimental and theoretical errors was obtained 
by Neubert [66] . This shows that there may not exist such a discrepancy. If we only employ 
the partial branching ratio for q 2 < 16 GeV 2 , a smaller \V u b\ will be obtained as given in 
Eq. ( 11021) . The analysis method is similar to the above and skipped here. 



D. q 2 shape 



BaBar has measured the differential rate for B — > nlui versus q 2 with good accuracy 
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531 ] . The q 2 shape of the differential rate relates to the q 2 shape of F^{q 2 ). It is 
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Figure 5: q 2 shape comparisons between C* form factor (the solid line) and experimental data. 
The dash line is a fitting to the experimental data. 

interesting to compare our prediction for the spectrum with the data. Based on F+ f , the q 2 
spectrum ABr/Br is calculated as shown in Fig. [5j We observe that the predicted shape of 
q 2 spectrum {C* q 2 shape) is inconsistent with the experimental q 2 spectrum. The \ 2 is 86.6 
for 12 degrees of freedom by assuming that the 12 bins of experimental data are completely 
uncorrelated. Parameterization methods [l4, 15], 1?1 18] have been widely used to determine 
the q 2 shape of F + (q 2 ) according to the q 2 shape data and theoretical inputs from LCSR and 
LQCD. It is an important task to compare F+ f (q 2 ) with these parameterization form factors 
to investigate their differences . This is left to other places. In summary, our calculations 
for the q 2 spectrum of the differential rate for B — > ttWi decays can not accommodate with 
the g 2 -spectrum of B — > nlui. The difference can be understood by the scaling behavior 
of F cf (q 2 ) and that of {dT/dq 2 ) th , Ff(q 2 ) oc 7T (1+aia0 and (T^dT/dq 2 )^ oc v ^ - 2 ^ for 
< q 2 < 10 GeV 2 and x = q 2 /{\GeV 2 ) 1 and Ff(q 2 ) oc r]~ 2 and (T^dT /dq 2 ) th oc r]- 1 for 
q 2 > 10 GeV 2 . We note that the form factor implied in BaBar q 2 data can have a scaling 
Fl A (q 2 ) oc r]^- 5 (x + r]y)^ 2 with x = -0.897 and y = 0.0281. This gives T^dT/dq 2 oc x + r] y . 
The residual scaling factor r] y in T~ 1 dT/dq 2 has small effects for large 77 (small q 2 ) and is 
important for small r] (large q 2 ). 
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VI. COMPARISONS WITH OTHER APPROACHES 



Many progresses in theories for B — > tx form factors were obtained in past years. Partial 



0(l/ma] corrections have been calculated in the hv factorization approach [3l|, |35] and 



LCSR te, 



T 

44| ] . The unquenched quark effects 10, 



11 



12 



13 



52] were included in the 



lattice QCD method. An all order proof of the factorization of the form factors was given 



for leading twist and twist-3 two parton contributions in the collinear factorization 



and in the soft collinear effective theory (SCET) [5 PL respectively. Different strategies for a 



49] 



parametrization of form factors have been developed 
our result with various theoretical methods. 



14 



15 



16 



17,08, 



3, 



231 ] . We compare 



A. Comparisons with k T 

It is the end-point divergent problem that the B — > it form factors contain logarithmic 
end-point divergences at twist-2 and linear end-point divergences at twist-3. To solve the 
divergent problem, the general wisdom is to use fc^-regularization for logarithmic end-point 
divergences. For linear end-point divergences, a threshold re-summation is needed, too. 
Twist-3 contributions dominate over the twist-2 ones in the kx factorization approach. This 
is known as a power expansion problem in k T . One possible solution to this problem is 
to employ twist-3 wave functions with better end-point behaviors. However, this would 
introduce additional uncertainties. 

For comparison, let's compare the values of F+ V,tw2 (0) calculated under the k R and £ R . 



Under k R , the result is F^ w2 (0) = 0.120 |3jj, where U = 0.13 GeV, f B = 0.19 GeV, and 
ujb = 0.40 GeV are used. Under £ R , F RlT ^ w2 (0) = 0.124 by using the same input parameters. 
This shows that £ R is as effective as k R at leading twist order. 

The large value -^+(0) is related to the power expansion method used in k T . In k T , only 
0(1/i71b) contributions from two parton Fock state of the pion were calculated. In |31| . 
<f)p{u) = 1 for the PS DA and <p a {u) = 6m (1 — u) for the PT DA of the pion were used. 



The collinear expansion scheme proposed by Beneke and Neubert (the BN scheme) 67 ] 
was employed. At q 2 = 0, the leading twist contributions lead to F_^ w2 (0) = 0.120 and the 
subleading twist contributions result in Ff. w3 (0) = 0.177. It is seen that F| w3 (0)/F| w2 (0) = 
1.5. In comparison, the £ R calculations give F^h(0) / F^? R (0) = 0.18. This shows that the 
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twist-3 contributions in fcj. dominate over the twist-2 contributions. On the other hand, the 



twist-3 contributions in C* are power suppressed as expected. The authors in [3lj proposed 
a new counting rule for the linear divergence associated with the twist-3 terms to explain 
this phenomenology. However, if higher order power corrections contain similar or even 
worse end point divergences, the power ex pan sion would breaks down. This is the power 



expansion problem in k T 35]. Hwang et al 35] pointed out that intrinsic parton transverse 
momenta are effective for the power expansion problem. The contributions related to the 
intrinsic transverse momenta are unknown, in principle. They proposed to use the model 

$(x,b)= d 2 k T e^{-ik T -b)^(x,k T ) (114) 

J\h\<l/b 

for all the wave functions (including twist-2 and twist-4 B meson WFs, and twist-2 and twist- 
3 two parton pion WFs) used in calculations, b are the conjugate coordinates to the intrinsic 
transverse momenta k^. Because the \l/(x,6) provides a better end-point suppression than 
the constant model for the twist-3 pseudoscalar pion WF, the twist-3 contributions become 
smaller than the leading twist ones, F+ w3 (0)/.Fj: w2 (0) = 1.5 — > 0.7. However, it is still unclear 
how higher order power contributions, such as twist-3 three parton and twist-4 contributions, 
can be included in the ktp approach. In addition, the introduction of intrinsic transverse 
momentum would arise a double counting problem when higher Fock state contributions 
are considered, such as three parton Fock state contributions, or four parton Fock state 
contributions. This is also unclear in kip. We argue that the power expansion problem 
in kp needs more efforts to clarify. Before this problem has been solved, the large value 
F7(0) = 0.297 in k T could be over-estimated. 

B. Comparisons with SCET 

In SCET, the form factor F+ n contain a factorizable and a nonfactorizable parts at 
leading twist order, where the factorizable part is expressed in terms of a nonperturbative 



form factor, 50] . By a fitting to the experimental data for the branching ratios of 
B — > 7T7T decays, (0) = 0.170 was found [36]. It is noted that our predicted value 

F + n (0) = 0.164 is close to the SCET result. This insures that the form factor Ff*(0) may 
not be as large as found from other approaches. The values derived by LSCR and LQCD 
and kp are about 0.26 — 0.30. The consistency between F+ f (0) and F^ CET (0) may not be an 
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accident. In the energetic limit, E 3> m n , the q(0)q(z) in the matrix element {0\q(0)q(z)\n) 
for defining the DA, becomes q(0)q(Xn / E) + ■ • ■ , where q(\n/E) is similar to the collinear 
quark field q c (Xn/E) defined in the SCET/. The collinear factorization based on full QCD 

m 

with energetic limit for light meson is likely equivalent to the SCET/ [49]. However, it needs 
further works to show the equivalence between C* and SCET/ for F^ n (q 2 ). 

The end-point divergences are also found in the SCET approach. The zero bin subtraction 
regularization method was proposed for dealing with soft and end-point divergences 3^]. 
Since the related physics about the subtracted quantities have not been given in the SCET 
language, there still remained some uncertainties in this method. 



C. Comparisons with LCSR 

The LCSR has been widely used for calculations of the form factors 



44. 
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691 ] . The LCSR method calculates f B F+ n (q 2 ) by matching the hadronic part to 
the partonic part of the correlator under parton-hadron duality. Let's first compare the 
LCSR form factor F^ CSR (0) and the form factor (0) at twist-2 order. Numerically, 



\F}: CSR ' tw2 (0)\ = 0.164 ~ 0.170 



44j and \Ff' tw2 (0)\ = 0.175 (lu b = 0.315 GeV) 



show that the LCSR and the C$ calculations are consistent at twist-2 order. At twist-3 
level, we find that the ratio of twist-3 contributions to the twist-2 contributions is equal to 
|F]: CSR ' to3 (0)|/|Fi CSR ' to2 (0)| ~ 80% - 100% fl, 3- Similarly, the ratio in the C> approach 
gives \Ff' tw3 (0)\/\F+ f ' tw2 (0)\ ~ 17% (u B = 0.315 GeV). We observe that twist-3 contri- 
butions in the LCSR form factor are more significant than those in the C* form factor. 
This also explain why the difference between \F° S ($)\ = \Ff' tw2 (0) + Ff' tw3 (0)\ = 0.164 
(w B = 0.315 GeV) and |FL CSR (0)| = ^^'^(O) + i^csR,^^ = a2? (at Q ^ ig 
so large. Note that there are partial cancellations in the sum of the twist-2 and twist-3 
contributions of the form factors. It was also noticed that the soft contributions from end 
point region in the LCSR analysis are important 0, [y]]. Because the twist-3 contributions 
have a more sensitive dependence on the end point behaviors of the pion DAs than the 
twist-2 ones, this may explain why the twist-3 contributions are so significant in the LCSR 
form factor. In LCSR calculation, the constant model for the PS DA of pion, <p p (u) = 1, is 
used. It is expected js, 31] that if the end point contributions can be properly dealt with 
by appropriate method, such as Sudakov factors, the soft contributions can be reduced. We 
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argue that the twist-3 contributions in the LCSR form factor could be overestimated. A 
better power expansion method may solve this problem. 



D. Comparisons with BN collinear expansion 



An expansion scheme developed by Beneke and Neubert (BN) 67| is widely used in lit- 
erature. The BN scheme is constructed for calculations of twist-3 two parton contributions 
from a final state pseudo-scalar meson. The BN scheme can not avoid end-point divergences 
in the hard spectator and annihilation contributions. There exists an ambiguity that the 
momentum and coordinate representations for the amplitudes are used in the calculations. 
Besides, the BN scheme is also used in the fey factorization 3l|, |35j . However, the distin- 



guishing differences between the collinear factorization and the kr factorization implies that 
these calculations require some cares. The reason is that the transverse parton momenta 
are assumed of order O(Aqcd) in the collinear factorization while they are not limited in 
the kx factorization. More detailed comparisons between the BN scheme and the collinear 
expansion refer to 41] . 



E. Comparisons with LQCD 

The lattice QCD approach can only calculate the form factors at large q 2 due to limits 
on the inverse space length of the ir meson energy. On the other hand, PQCD approach is 
applicable for small q 2 , where the virtual radiations are perturbative. Due to £ R , the 
approach developed in this paper for the form factors has extended the application range 
from low q 2 to moderate q 2 . These two approaches are complementary and can be combined 
to derive form factors of full range q 2 . The related works are left in our future publications. 
The comparisons between the extrapolation of F? (q 2 ) to large q 2 region with the LQCD 



calculations by FNAL 



11 



12 



52] and HPQCD [13] collaborations show that the C* 
prediction is close to the FNAL calculations (with a x 2 /N = 54.8/15) but has a large 
deviation with the HPQCD calculations (with a x 2 /N = 271.1/7). The C* prediction is 
close to the FNAL calculations than the HPQCD calculations. 



40 



VII. DISCUSSION AND CONCLUSIONS 



In this paper, we have proposed a ^-regular ization for the logarithmic end-point diver- 
gences in B —>■ 7i form factors. It can effectively resolve the end-point divergent problem. 
The linear end-point divergences are solved by an simultaneous use of the ^-regular ization 
and non-constant twist-3 distribution amplitudes. The £-regularization and the existence 
of a non-constant pseudos-scalar DA are followed by the collinear expansion. The factor- 
izability of B — > it form factors is shown up-to 0(ot s /rriB) by including complete twist-3 
contributions from the 7r meson and partial 0(l/m 2 B ) from the B meson. The calcula- 
tions are given explicitly. The form factors are calculated and applied to determine the 
CKM matrix element \V u b\ according to the world average branching fraction. Eq. fllOip 
is our main result. The result is also applied to make a prediction for the q 2 spectrum of 
B — > ixlvi. The predicted differential rate is inconsistent with the BaBar spectrum data. 
This discrepancy deserves further studies. Generalization of our result to other form factors 
are straightforward. We leave these interesting tasks to other places. 

The twist-3 contributions calculated in the h T and LCSR approaches have been compared 
to our calculations. We argue that the twist-3 contributions in these two methods could 
be overestimated. If the twist-3 contributions in these two methods can be reduced by 
an appropriate power expansion method, it is possible to have a consistent result with 
our calculations. However, before the method has been derived under the k T or LCSR 
formalisms, the k T and LCSR form factors are larger than the form factor. 

One important application of the present work is to improve the QCD factorization at 
leading order in l/m& expansion. The QCD factorization at leading order in expansion 
demonstrates that the matrix element (P\P2\Qi\B) for B — ► PP decays can be expressed 
by the factorization formula 



(P 1 P 2 \Q i \B) = F BP > / di / du^Tl (£,ui)0 B (£)MM] + (Pi <- P2) (H5) 



where the T I,n are the type-I and type-II hard scattering kernels for non-factorizable radia- 
tive corrections. Due to the end-point divergences in the form factor F BPl and in twist-3 
contributions, the factorization formula is only valid at leading twist (twist-2). The form 
factor F BPl is isolated from the factorization formula and identified to be a physical form 



+ / d£ / dm / du 2 Tr[T i (f, u u M2)0b(O0p 2 («i)0p 2 («2)] , 
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factor, which can only be determined by experiments. According to the results derived in 



this paper and those obtained in 41(, we propose a generalized QCD factorization formula, 
which is valid at twist-3 order. The generalized QCD factorization formula for B — > PP 
decays with PP two light pseudo-scalar mesons is given (under two parton approximation) 

(PMQilB) = f P2 j dij duxTrp^^MO^M] + (Pi <- P2) (116) 
+ / "df / "dm j du 2 TT[T l I ^,u 1 ,u 2 )(pB(04>Pi(u 1 )4>p 2 (u2)} 

+ [ ^ f dUl j du 2 Tl [ T i I (^ U l^ U 2)(pB(0 ( t>P2( U l) ( t>P2( U 2)], 

where T° is the hard scattering kernel for factorizable radiative corrections, T 1 is the hard 
scattering kernel for type-I non-factorizable radiative corrections, and T n is hard scattering 
kernel for type-II non-factorizable radiative corrections. The £-regularization is used to 
regularize the end-point divergences in form factors, the twist-3 annihilation contributions, 
and twist-3 hard spectator contributions in charmless B — > PP decays 4JJ. The leading 
order is 0(a s ). The improvement in this generalized QCD factorization is that the form 
factors are calculable and the precision order is generalized to include 0(1 /ms) corrections 
under two parton approximation. The collinear expansion is a necessary tool to derive 
correct twist-3 contributions in the above factorization formula. The application and the 
formal proof of the generalized QCD factorization formula will be given elsewhere. 

The generalized QCD factorization formula is valid at twist-3 under two parton approxi- 
mation. The factorization formula of the same twist order have also been shown to exist in 
the /ct factorization and in the SCET, respectively. A consistent picture from three different 
approaches is obtained that the factorization theorem for B — > PP decays is valid up-to 
twist-3 order under two parton approximation. 

Appendix A: REDUCED HARD SCATTERING FUNCTIONS 

The reduced hard scattering function H$ (£, u) is derived according to the rules given in 
Section II. According to the Feynman diagrams depicted in Fig. 1 (a) and (b), the reduced 
hard scattering function H^\^,u) is H^\^,u) = H^'^ a \^,u) + Hp'^ b \£, u) where 
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C b 7ra s 



(Al) 



W ),(6) (e, u)] ij;kl = (2eM, i [^] w -^[^ 7a]i .[ 7 ,]^f)^. 

(A2) 

The reduced hard scattering function Hj£f v (£ ,u q ,u q , u g ) are derived according to the Feyn- 
man diagrams depicted in Fig. [3)^1) - (4). The contributions from the other diagrams de- 
picted in Fig. EI5) - (36) are of higher twist order than twist-3. They are neglected in 
this work. The expression for Hj$ v (£ ,u q ,u q ,u g ) is also written as H$ v (£ ,u q ,u q ,u g ) = 
H^ {a) u q , u g ) + H^ b) (t,u qi u q , u g ) where 

(A3) 

[H^Ln q ^u g )l m =■ -°* ^ 



N c plr] 2 E A u q u q u g ^ - 7]u q ) 



x ( 2£ [lul v fiL [i\ k i- r)U q [ 7 „ 7?) i fi 7a ] [ 7/9 ] w df ) p W)A , . (A4) 



Appendix B: 5 MESON DISTRIBUTION AMPLITUDES 



When sub leading order effects in l/m# expansion are considered, the B meson DAs 
^b" 2 (0 an d 0b (0 are defined as [46j 



(01^(^(0)15) = !^(( y p B + ms ) / 



where 0s' 2(7 ' ) (O and 0b (0 correspond to the 0#(£) and t(£) — 0b (0 defined in the liter 



4iV c 



>g*M(£) +£(,(_ 7s)(Bl) 



ature. The models for 0j(£) and B (£) are assumed to be j46j 



0+(O = ^?£exp 



B 



B 



(0 



^B 



L 



exp( £) 

UJ B oj b 



Their first two inverse moments should satisfy the conditions 



fd^'-\0 = 1 
Jo 



i A+,-), 



,^ 0b' (0 ^ < 3 mg 



(B2) 
(B3) 

(B4) 
(B5) 
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where A = mg — m&. In this work, we employ the following models for </>^ 2 ^ and 4>b 

0b 2(/) (O = <t>UO = ^eM-—0, (B6) 
MO = ««) - = f ({ - «P<-^« - ^> , (BT, 

where 9 is a step function to insure that 0_b(£) is positive. 
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